UNCLASSIFIED 


405  720 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 

Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


Z  l<ZOf 


FTD-TT* 


63-28^ 


405  720 


63-3- 


TRANSLATION 


turbulent  boundary  layer  ik  a  gas  flow  with  heat  transfer 
and  a  prakdtl  number  different  from  unity 


Yu.  V.  Lapin 


FOREIGN  TECHNOLOGY 
DIVISION 


AIR  FORCE  SYSTEMS  COMMAND 


WRIGHT-PATTERSON  AIR  FORCE  BASE 
OHIO 


■mM 


\h  i~,  t-v 

TTinnn  np 

•M.'1  ■  10<jg 


UNEDITED  ROUGH  DRAFT  TRANSLATION 


TURBULENT  BOUNDARY  LAYER  IN  A  GAS  FLOW  WITH  HEAT  TRANSFER 
AND  A  PRANDTL  NUMBER  DIFFERENT  FROM  UNITY 

BY:  Yu.  V.  Lapis 

English  Pages:  15 

SOURCE:  Russian  Book,  Tekhnicheskaya  Gldromekhanlka,  Gosudar- 
stvennoye  Nauchno-Tekhnicheskoye  Izdatel'stvo  Maahln- 
03troitel'noy  Llteratury,  Trudy  LPK,  No.  217,  Moskva, 
1961,  pp.  27-36 


S/563-61-0-217 


THIS  TRANSLATION  IS  A  RENDITION  OR  TNE  OR»t> 

NAL  P0REI6N  TEXT  WITHOUT  ANY  ANALYTICAL  OR 
EDITORIAL  COMMENT.  STATEMENTS  CR  THEORIES 

ADVOCATED  OR  IMPLIED  ARE  THOSE  OR  THE  SOURCE 

AND  DO  NOT  NECESSARILY  REELECT  THE  POSITION 

TRANSLATION  DIVISION 

OR  OPINION  OP  THE  POREION  TECHNOLOGY  OL 

POREIGN  TECHNOLOGY  DIVISION 

VISION. 

WP-APR.  OHIO. 

FTD-TT-  63-284/1+2+4 


p 


•  TURBULENT  BOUNDARY  LAYER  IN  A  GAS  PLOW  WITH  HEAT 
TRANSFER  AND  A  PRANDTL  NUMBER  DIFFERENT  FROM  UNITY 

Yu.  V.  Lapin 

The  determination  of  surface  friction  and  of  heat  transfer  in  a 
turbulent  boundary  layer  during  the  movment  of  substances  with  large 
supersonic  velocities  is  an  urgent  problem  of  aerodynamics.  By  now 
there  have  appeared  a  significant  amount  of  articles  devoted  to  this 
problem.  However,  the  more  or  less  full  consideration  of  the  influenci 
of  all  factors  (pressure  gradients,  radiation,  and  others)  on  the 
dynamics  and  heat  characteristics  of  a  boundary  layer  remains  very 
intricate. 

The  existing  methods  of  calculation  of  turbulent  boundary  layer 
in  a  gas  with  heat  transfer  for  both  the  case  of  a  flow  past  a  flat 
plate  and  for  the  flow  with  pressure  gradients  are  based  on  the  trans¬ 
fer  to  gas  dynamics  of  the  formulas  of  semiempirical  theories  of  turbu¬ 
lence.  This  article  is  In  this  sense  no  exception.  The  solution  is 
obtained  for  the  flow  of  gas  with  a  moderate  pressure  gradient  and 
with  an  arbitrary  temperature  distribution  at  the  wall  since  the 
Kannan  semiempirical  theory  is  used.  The  influence  of  the  pressure 
gradient  is  calculated  only  across  the  external  flow.  The  temperature 
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range  Is  examined  in  which  the  Prandtl  number  and  the  specific  heat 
can  be  considered  constant,  regardless  of  temperature. 

In  this  article  we  obtain  the  analytical  expression  for  the 
parameter  H  »  &*/£**  depending  on  the  conditions  at  the  wall  and  at 
the  exterior  boundary  layer.  The  final  formula  obtained  in  this 
article  for  the  calculation  of  friction  does  not  require  the  use  of 
various  types  of  auxiliary  graphs  and  tables,  as  in  other  existing 
methods.  This  advantage  can  prove  to  be  beneficial  in  a  number  of 
cases.  The  calculation  of  heat  transfer  is  based  on  Reynold's  analogy. 

1.  Basic  equations 


The  differential  equations  of  continuity,  of  momentum  and  of 
average  stationary  plane  kinetic  energy  in  a  turbulent  boundary  layer 
have  the  form: 


du  .  du 

~3y  ~ 


dp  d  \ .  du  1 

J: 
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(1) 

(2) 
(3) 


where  x,  y  and  u,  v  —  are  respectively  the  coordinates  and  the  compo¬ 
nents  of  the  velocity  along  the  surface  and  normal  to  it; 

p  *  density  of  the  gas; 
p  *  pressure; 

p  and  e  are  respectively  the  coefficients  of  dynamic  and  of  eddy 
■  ~  viscosity; 

T0  *  stagnation  temperature  in  a  boundary  layer; 

T  *  temperature  of  the  gas  in  a  boundary  layer;  ,  -  -  — 

’  Cp  *  specific  heat  with  constant  pressure;  1 

:-r;-  *  Yrc y  hi*? 
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J  *  the  mechanical  equivalent  of  heat; 

f  '•  ■  ■ 

g  -  the  accelerating  foi*ce  of  gravity; 

pram-Bjt.  m  the  Prandtl  number; 

X  **  the  coefficient  of  heat  transfer  of  the  gas. 

The  boundary  conditions  are: 

u—v—  0,  T  —  Tw  (a)  when  y  —  0;  1 

Ut  (x).  T  =  T,  (x)  when  y  =  co.  )  ' 

We  will  consider  the  external  flow  as  lsoentropic,  i.e.  the  total 
pressure  pgo  and  stagnation  temperature  Tgo  are  constant. 

Outside  of  the  boundary  layer  Bernoulli's  equation  holds  true: 


—0U  dU' 

.  dx  fa  • 


(5) 


Prom  equations  ( l)  ,  (2),  and  (5)  and  boundary  conditions  (4) 
the  integral  relation  of  momentum  can  be  obtained  by  the  usual  method: 


— C-  P(x)6**  =- 1  T‘ 

dx 


A*  T~ 


Here 


where 


P(x)  = 


2  +  H  —  Ml  dM . 


H 


b% 


* 


=  J  — Is  the  thickness  of  the 
0  * 

momentum  loss; 

m 

6*=J‘^l — j  dy  _  is  the  displacement 
0 

thickness, 

h  -  «  t/~^T 


(6) 

(7) 

(8) 

(9) 

(10) 

(11) 
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where 


r(p  =  pa  —  is  the  friction  at  the  wall; 

c.  =  —  is  the  coefficient  of- -fric-tlon. 

o  y,  ■  • 

Here  and  throughout  the  article,  a  dash  will  designate  the  dimen¬ 
sionless  values  (  x  ®  =  — t~  >  u  =-w-  ,...),  obtained  by 

ue 

division  of  the  dimensional  value  by  a  selected  scale  (for  example, 

L  -  airfoil  chord) . 

2.  Relation  Between  Velocity  Profile  and  Temperature  Profile 


In  the  laminar  sublayer,  the  expression  for  stagnation  tempera¬ 
ture  will  be  sought  as  a  series  in  powers  of  uj  in  view  of  the  small¬ 
ness  of  u  ( separation  close  to  the  wall)  we  will  confine  ourselves 


to  three  terms 


T#  •*  fio  +  axu  +  CjU*. 


We  will  note  that  To  *  T^  when  u  =  0  and,  consequently, 

a.  =  T„ 


(12a) 


Having  differentiated  equality  ( 12)  with  respect  to  jr  and  having 


taken  the  obtained  expression  when  y  -  0,  we  will  find  that 


(12b) 


For  the  determination  of  the  coefficient  a2  we  will  twice 
differentiate  equality  ( 12)  with  respect  to  £  and  we  will  use 
equations  (2)  and  (3)  where  y  “  0;  then  after  simple  conversion. we 


will  obtain 


1  —  Pr  _  Prq~uw  dp 

2SJc„  dx 


yc?  ,-f-£ 


(_12c) 


4. 


By  the  examined  mean  pressure  gradients,  the  second  term  in 
equality  ( 12)  proves  to  be  negligible  in  value  and  can  be  omitted  in 
the  future. 

We  will  introduce  with  respect  to  h  the  value  h. ,  expressed  by 
the  equality  r; 


(13) 


where  the  Stanton  number 


S/  _ _ v* 

QeUtgCp(Tt  —  Tv) 


(it) 


represents  the  dimensionless  coefficient  of  heat  transfer,  *=  —  (*■■§£■). 
is  the  heat  flow,  is  the  equilibrium  temperature  of  the  flow  sur¬ 
face  without  heat  transfer. 

Using  equality  (ll)  and  the  expression  for  stagnation  teraperatur* 


(15) 


we  will  convert  the  relation  ( 12)  to: 

r  _ 

~7T  -  1  ~  Prau — Pryu2, 


(16) 


whei*e 


(17) 


In  the  turbulent  core  we  will  use  Reynolds  analogy,  which  in 
the  case  of  a  compressible  gas  is  expressed  by  the  formula  [l]: 


+  <l8> 

If  we  simply  assume  the  constancy  of  the  friction  force  x  * 
and  oi‘  heat  flow  q  =  q^  across  the  boundary  layer  and  integrating 
( 18)  we  will  have:.  . “  '  -  V..,:  - 


(19) 


T°  u  -rC  (a). 

Determining  C(x)  from  the  condition  at  the  boundary  of  the  laminar 
sublayer  we  will  obtain  after  conversion  the  following  relation  of 
temperature  profile  with  velocity  profile  in  a  turbulent  core: 


~t~  —  l  —  P  —  ««  —  VM\ 

*  w 


(20) 


where 


P  =  (1  —  Pr)  (o)Ujt  --  y u\^. 


(21) 


and  ult  is  the  dimensionless  velocity  at  the  boundary  of  the  laminar 
heat  sublayer. 

Prom  equality  (20)  and  boundary  condition  T  =  Te  where  u  =  1 
we  will  obtain  the  expression  for  co: 


(I)  : 


i - 


(22) 


The  previously  obtained  relations  between  temperature  profile  in 
a  laminar  sublayer  (  16)  and  in  a  turbulent  core  ( 20)  in  the  case'  of  a 
flat  plate  =  0)  and  Pr  =  1  convert  to  the  known  Crocco  relation: 

(23) 

Prom  the  equalities  ( 16)  ,  (20)  and  the  equation  of  the  component 


p  =  qRT 


(24) 


the  expression  for  surface  gas  in  a  laminar  sublayer  is  obtained: 


o  _ . 

Ci 


1  -f  ft  —  a  —  Y 


1  —  Pr  —  Pryu 


(25) 


and  in  a  turbulent  core: 
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\ir 


(26) 


1  ~  ft  —  n>  —  y  ^ 

|  —  —  tuxv  —  yu* 


J>.  Determination  of  Integral  Thicknesses  of  a  Boundary  Layer 

We  will  show  the  method  of  determining  the  Integral  thicknesses 
by  the  example  of  the  calculation  of  momentum  loss  thickness  6**. 

In  the  expression  for  momentum  loss  thickness  ( 19)  we  will  con¬ 
vert  to  the  universal  coordinates: 


then  we  will  have: 


6**  =  — f  JL .  V  ( ,  <r  \  , 

We  can  present  equality  (28)  also  In  the  form: 


*-  J-J-SO-BW. 


(27) 


(28) 


(29) 


where 


The  relations  of  densities  -  in  a  laminar  sublayer  and  in 
a  turbulent  core  are  expressed  by  equalities  (25)  and  (26)  respect¬ 
ively.  Thus,  for  the  determination  of  6**  it  is  necessary  to  have  the 
relationship  of  the  value  i)  to  dimensionless  velocity  u. 

In  a  turbulent  core  for  the  determination  of  tj  we  will  use  the 
Karman  relation! 

T  =  CxsfLj,  (ZO) 

_  :  U  ~ 

where  x  is  the  constant  of  turbulence,  the  value  of  which  is  taken 
equal  to  x  -  0.4  as  in  the  case  of  an  Incompressible  fluid;  the  prime 
in  equality  (30)  designates  a  derivative  with  respect  to 

Converting  in  (30)  to  universal  coordinates  and  bearing  in  mind 

_  .  -7r 


equality  (20)  after  simple  conversions  we  will  obtain  the  equation 


V  i  T  p  —  < 


yu- 


(31) 


The  minus  sign  appears  after  the  root  extraction  because  V<  0, 

If  the  Influence  of  the  pressure  gradient  on  the  velocity  profile 
is  not  considered. 

In  equality  (}l)  let  us  interchange  the  argument  and  the  function, 
then  we  will  obtain  the  equation: 

j\  _ _  * 

}  1  ~  p  —  oiu  —  y5*  (32) 


the  equality 


t)  =  C  exp 


v'y  a  ■ 


2Vy 


V 


(33) 


is  used  for  the  first  integral  of  (32)  . 

For  the  determination  of  the  constant  C,  we  assume  that  this 
value  is  the  same  as  in  the  case  of  incompressible  gas  on  the  flat 
plate,  i.e.  C  =  -jr  *  wa,  where  a  is  the  turbulent  constant,  the  value 
of  which,  as  in  the  case  of  incompressible  gas,  is  taken  to  equal 
a  =  11.5.  In  other  words,  it  is  assumed  that  the  value  of  the  deriva¬ 
tive  q  in  the  case  of  a  turbulent  core  does  not  depend  on  compressa- 
bility  and  gradient  pressure.  We  note  that  the  assumption  of  the 
non-dependence  of  the  value  tj  on  compressability  with  no  pressure 
gradient  agrees  with  the  experiment. 

The  assumption  leads  to  the  following  result! 


q  =  ~r  exp 


Vy  “+ 


arc  sin 


2Vy 


*  T  P  + 


-  arcsu 


_  2  Vy 


V  »+p 


+ 


(34)  — 


where  o  *  ^  ^ 
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The  Integral  in  the  right-hand  part  of  equality  (29),  strictly 
speaking,  should  be  divided  into  two  parts:  0  <  u  <  a  and  a  <  u  <  1 
and  could  substitute  their  own  value  for  the  values  n  and 

IV 

In  each  of  these.  However,  in  view  of  the  relative  thinness  of  the 
laminar  sublayer,  we  can  omit  the  first  part  having  continued  the 
second  up  to  the  wall. 

Thus,  excluding  tj  with  the  help  of  equality  (34)  and  substitutinj 
the  value  o-'o,  from  ( 26)  ,  we  will  have  the  following  expression 
for  the  momentum  loss  thickness: 

i 

a  *  *  _  j|  P — uj — y)  i*  u  n  *j  _  ^ 
lU.  J  i  j-  p  —  u u — yu* 


X  exp 


—z  I  arc  i>m  — 7=============-=- 

Kv‘  ]/i-H>+-£ 


Vy« 


—  arc  sm  ■ 


2  V  Y 


V‘+»+-jr 


du. 


Introducing  a  new  variable 


V\i+ 


we  will  obtain 


V  =  r^arc  sin  — 27(X— . 


(l  +p_tt__v)VcA!  1/ J  +I|  + J£-  f  . 

1  **  - - — ! _ iv_  1  e*p[y~i<v—'pJi  v 

^  .)  cas 


where 


X  [sin  ^  -  sin  ) fy  xj-J  [sin  )  ^  t,h,  —  sin  ]/y$]  dxp. 

Kv-}-  — 


“ipg  =  17=-  arc  sin 
Vr 


2  V Y 


K'+t+4-  • 


I 

-7^-  arc  sin 
n’ 

2  t^y 

1 

*-r=r  arc  sin  - 

—  —  .  (a 

1'  v  Cl  -A 

"  1  u 

_  2J/Y 

kY 

V‘+f+-W 

;t, P 
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(35) 


(36) 

(37) 


(38) 


The  value  Tih,  as  Is  obvious  from  Its  determination,  essentially 
exceeds  unity,  therefore  for  the  calculation  of  Integral  (37)  we  use 
a  general  representation  of  the  Integral  containing  an  exponential 
function  In  the  form  of  an  asymptotic  series  which  is  obtained  as  a 
result  of  integration  by  parts: 


•  r  1  1  1 

j  .~'-v  f  (v)^ = [/  (V) — -ar  /*  (*>  +  r  W) — •] 


(39) 


Fulfilling  in  (37)  the  integration  accurately  to  the  term 
containing  x3h3  in  the  denominator  and  representing  approximately 


as 


m 


we  will  obtain  the  following  expression  for  the  momentum  loss 
thickness: 


6**  = 


*~*°  vtt  (i  j-  p — —  y)  „** (*«-*») 
T^u. 


/.  2-^-l,5a  — y  \ 

\  *kY~i  (4l) 


The  following  expression  for  the  displacement  thickness  can  be 
analogously  obtained: 


6*  -  *  X°v--  d-f  P-fV)  x 

jVrU, 


Xfl 


y(3  4-30-1,50  — 


where 


*AO-rP-i-Y)Kl+P-“-Y 


Vy  -r 


2Vy 


ib,  —  =  —r=  arc  sin - > - 

l/.+HH-g- 


-arcsm  ■ 


2Vv 


]/  H-P-f 


— 

*Y 


(42) 


(43) 


If  we  set  up  the  ratio  of  the  equalities  (42)  and  (4l)  we  will 
obtain  the  expression  for  H: 
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•-  f  V  -r  3j5  —  l,5w  —  y  -f.  \  1 

(1 T  P  T  V)  1 - - - - - ( 44) 

//  _ I _ x/,  (1  —  jt  --  yi  Y  1  +•  p  —  in  —  y  J  '  ' 

(1  —  P  —  <d  —  V«  1 - ~r‘— -  - 

\  *AJ  1-H»  — to  — Y  / 

In  a  limiting  case  of  the  flow  of  an  incompressible  gas  without 
heat  transfer  we  will  obtain  the  following  expression  for  H: 

■■  "“Tibr*  (*5) 

yJi 

which  with  Reynolds  numbers  10®  4-  107  yields  the  value  H»1.3>  Which 
agrees  with  the  experiment. 

4.  Solution  of  Integral  Equation  of  Komenta 


Introducing  a  new  variable  in  the  equation  of  momenta  (6) 


=  3  *  *  exp  f 


P  (x)dx  \  *. 


we  will  reduce  It  to 


where 


dZ  _  T, 
<u 


N  ~  f  P(x)dx, 


and  Is  the  dimensionless  coordinate  of  the  beginning  of  the 
turbulent  part  of  the  boundary  layer. 

Integrating  equation  ( 47)  and  reverting  6**  we  will  obtain: 


=*■"()  4r~frdx+c 


If  in  (49)  in  place  of  7**  we  substitute  its  value  from  equality 

(4l)  we  will  have:  „  — — 

2  +  2P-l^-V  _ vl  f  /_  .  —  - 

f(A)l  *kYl+t-m^r  ~  **  T*>  J.  (50) _ 
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where 


Plrt  _  j  xg(l  -ip  — a>  — v)vg 
w  ~  i*'UA- 


(51) 


Equation  (50)  permits  us  to  use  the  method  of  successive 
approximation  for  the  determination  of  h.  Carrying  out  the  conversion, 
we  will  obtain: 

/■£-•-£■*+«  (52) 

7  1  ,  _ % _ 

1  V.  dr  —  *«.>  /j 1_\  ‘ 


The  constant  of  integration  C  is  determined,  as  usual,  from  the 
condition  of  the  congruence  of  momentum  loss  thicknesses  at  the  point 
of  conversion  of  the  laminar  boundary  layer  to  the  turbulent  boundaxy 
layer.  Upon  completion  of  the  turbulent  boundary  layer,  the  constant 
of  integration  C  is  reduced  to  zero,  and  equality  (52)  is  converted 


to 


h  = 


_ I _ 

MV*— 


In 


(I 


£  -  -ft  —  —  Y  I  ’ 

y.h  J'l-p-w-yj 


(53) 


It  Is  natural  to  take  the  value  h  calculated  for  a  flat  plate 
as  a  zero  approximation. 

5.  Calculation  of  Heat  Transfer 

Using  the  equalities  ( 17)  and  ( 22) ,  we  will  find  the  expression 

for  the  parameter  (-^— )  *s 

t 

/_M*_  2St  l-yf  +  O-Pr)  v?x±  (54) 
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The  numerator  of  the  right-hand  part  of  equality  (5^)  can  be 

Tt 

considered  equal  to  1  —  •**-  although,  strictly  speaking,  it  differs 
somewhat  from  that  remainder  since,  in  the  conversion  to  the  case  of 
^o)  *  the  value  u^t  changes?  however,  the  influence  of  this  change 
is  negligible.  Therefore,  with  great  accuracy  the  equality  (5*0 
can  be  written  as: 


2St 


(55) 


cl  1  —  (1  —  Pr)\% 

This  formula  agrees  in  form  with  the  Prandtl  formula  [2], 
obtained  by  him  for  the  flow  of  an  incompressible  fluid  in  a  tube. 

The  value  of  a  dimensionless  velocity  at  the  boundary  of  a  heat 
laminar  sublayer  is  expressed  by  the  equality: 


(56) 


where  °t  =  :t=-  is  the  constant  of  a  heat  laminar  sublayer. 

V  Pr 

We  note  that  formula  (55)  yields  the  results  agreeing  with  the 
A.  R.  Colburn  experimental  data  \3\>  who  found  that  with  small 
velocities  and  minor  temperature  differences  the  relation  of  the  Stanton 
number  to  one-half  of  the  friction  is  expressed  by  the  ratio: 


6. 


Calculation  Example. 


(57) 


Calculations  of  friction  and  heat  transfer  on  a  wing  profile  in 
a  supersonic  flow  were  carried  out  by  the  stated  method.  The  profile 
described  by  the  formula  y  =  0.l6  3c(l  —  x)>  at  angle  of  attack  0  =  5° 
was  taken  as  the  initial  profile. 

The  distribution  of  numbers  on  the  profile  were  calculated  — 
with  respect  to  the  conventional  formula  of  the  theory  of  a  profile  in 
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supersonic  .flow. 


1  ~o.nji 
j/Aii-i 


8  + 
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0,0-Ull  —  —  0.4.M*  —  0.5  •  . 

- - 7^ - V  . 

where  0_  is  the  local  angle  of  attack. 

P 

In  the  external  flow,  conditions  at  a  height  of  flight  H  =  20  km 
where  6  were  used.  The  temperature  at  the  wall  was  taken  to  be 

Tw  -  700°  K.  . 

The  results  of  the  calculation  of  the  coefficient  of  friction 
c^  and  of  the  parameter  2St/c^.  are  presented  in  the  diagram.  As 
appears  in  the  diagram,  for  practical  purposes  we  can  be  satisfied 
with  the  results  of  the  second  approximation. 
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